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In this article, we present an emerging field of quantum chemistry 
at finite temperature. We discuss its recent developments on both 
CN . theoretical and experimental fronts. We describe and analyze several 

experimental investigations related to the temperature effects on the 
5o ' structure, electronic spectra, or bond rupture forces for molecules. 

This includes the study of the temperature impact on the pathway 
shifts for the protein unfolding by atomic force microscopy (AFM), 
Tj- ' the temperature dependence of the absorption spectra of electrons 

■ in solvents, and temperature influence over the intermolecular forces 

^ , measured by the AFM. On the theoretical side, we review a recent 

• 1— I ' advancement made by the author in the coming fields of quantum 

chemistry at finite temperature. Starting from Block equation, we 
, have derived the sets of hierarchy equations for the reduced density 
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operators in both canonical and grand canonical ensembles. They 
provide a law according to which the reduced density operators vary 
^ . in temperature for the identical and interacting many-body particles. 

^ I By taking the independent particle approximation, we have solved the 

equation in the case of a grand canonical ensemble, and obtained an 
eigenequation for the molecular orbitals at finite temperature. The 
explicit expression for the temperature-dependent Fock operator is 
also given. They will form a foundation for the study of the molecular 
electronic structures and their interplay with the finite temperature. 
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Furthermore, we clarify the physics concerning the temperature ef- 
fect on the electronic structure or processes of the molecules which is 
crucial for both theoretical understanding and computational study. 
Finally, we summarize our discussion and point out the theoretical 
and computational issues for the future explorations in the fields of 
quantum chemistry at finite temperature. 

Keywords Quantum chemistry at finite temperature; temperature depen- 
dent; polymers; protein folding; intermolecular forces; solved electrons. 

1 Introduction 

The history for quantum chemistry development is almost synchronous to 
that of quantum mechanics itself. It begins with Heitler and London's study 
of electronic structure of H2 molecule shortly after the establishment of wave 
mechanics for quantum particles There are two major types of molecu- 
lar electronic theories: valence bond approach vs. molecular orbital method 
with the latter being the popular one for the present investigation. It has 
gone through the stages from the evaluation of molecular integrals via a 
semiempirical way to the one by an ab initio method. Correlation issue is 
always a bottleneck for the computational quantum chemistry and is un- 
der intensive study for over fifty years. For large molecular systems such 
as biomolecules and molecular materials, the development of the combined 
QM/MM approach, pseudopotential method and linear scaling algorithm 
has significantly advanced our understanding of their structure and dynam- 
ics. There are about eight Nobel prize laureates whose researches are related 
to the molecular electronic structure theory. This not only recognizes the 



2 



most eminent scientists who have made the outstanding contributions to the 
fields of quantum chemistry, but more importantly, it indicates the essen- 
tial roles the electronic structure theory has been playing in the theoretical 
chemistry as well as for the whole areas of molecular sciences. Nowadays, 
quantum chemistry has been becoming a maturing science [21 E] ■ 

Nevertheless, the current fields of quantum chemistry are only part of the 
story for the molecular electronic structure theory. From the pedagogical 
points of view, the quantum mechanics based on which the traditional quan- 
tum chemistry is built is a special case of more general quantum statistical 
mechanics jU IH] . In reality, the experimental observations are performed 
under the conditions with thermodynamic constraints. Henceforth, there is 
a need to extend the current areas of quantum chemistry to the realm of, for 
instance, finite temperature [11131^1. 

Indeed, many experimental investigations for various fields and for dif- 
ferent systems have already shown the temperature or pressure effects on 
their microscopic structures [7-30,49-57,63,71-87]. The polymeric molecules 
are one of the most interesting systems for this sort of studies [7-16]. The ex- 
perimental measurement on the absorption spectra, photoluminescence [PL), 
and photoluminescence excitation {PLE), and spectral line narrowing {SLN) 
for the PPV and its derivatives all show the same trend of the blue shift 
with the increasing temperature [3 El El- This attributes to the temperature 
dependence of their very rich intrinsic structures such as the vibronic cou- 
pling IT^l IT^ . The experimental investigation of the temperature effect 
on the biomolecules started in the late nineteenth century [T71 118j. Most re- 
cently, it has been extended to the study of folding and unfolding of protein 
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or DNA |2ni 1^ • In addition to the observed patterns for the unfolding 
forces with respect to the extension or temperature, it has been proved that 
the temperature-induced unfolding is another way for the study of mech- 
anisms or pathways of protein folding or unfolding processes [19-23]. The 
newest related development is on the AFM measurement made by Lo et al. 
of the intermolecular forces for the biotin-avidin system in the temperature 
range from 286 to 310K [US]. It has shown that an increase of temperature 
will almost linearly decrease the strength of the bond rupture force for the 
individual biotin-avidin pair. The study of temperature effect on the absorp- 
tion spectra of solvated electron began in the 1950 's and it is still of current 
interest. A striking effect is that an increasing temperature will cause the 
positions of their maximal absorption red shift [71-85]. 

In recent papers, we have deduced an eigenequation for the molecular 
orbitals [HEI- It is the extension from the usual Hartree — Fock equation at 
zero temperature to the one at any finite temperature jHHllHni- It opens an 
avenue for the study of the temperature effects on the electronic structures 
as well as their interplay with the thermodynamic properties. In the third 
section, we will present this equation and give the details for its derivation. 
In the next section, we will show four major types of experiments related 
to the study of the temperature influences over the microscopic structure 
of molecular systems. In the final section, we will discuss and analyze our 
presentations, and point out both theoretical and computational issues for 
the future investigation. 
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2 Experimental Development 



In this section, we mainly describe the experimental investigations related 
to the temperature effect on the bonding, structure and electronic spectra of 
molecules. We choose four kinds of the most recent developments in these 
fields which are of chemical or biomolecular interests. 

2.1 Temperature effects on geometric structure and 
UV-visible electronic spectra of polymers 

The first important systems where the important issues related to the tem- 
perature effect on the geometric structure and electronic spectra are the 
polymeric molecules. Many experimental investigations and some theoreti- 
cal work already exist in the literature [7-16]. However, how the temperature 
changes the microscopic structures of the polymers are still not completely 
understood and there are many unresolved issues in interpreting their elec- 
tronic spectra. We list here a few very interesting experimental investigations 
for the purpose of demonstration. 

The poly(;5-phenylenevinylene)(PPF) is one of the prototype polymeric 
systems for the study of their various mechanical, electronic, and optical 
properties. The impact from the temperature on the absorption spectra, 
photoluminescence (PL), and photoluminescence excitation [PLE) of the 
PPV have also been investigated both experimentally and theoretically [71 
IHlini- In an experiment carried out by Yu et al., the absorption spectra are 
measured for the PPV sample from the temperature 10 to 330-ft". The details 
of the experiment are given in their paper jH] . The resulting spectra for the 
absorption at T = 80 and 300K are shown in Figure 1 of that paper. We 
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see that there is a pronounced change in the spectra when increasing the 
temperature. They also study the PL and PLE spectra for the PPV. The 
measured PL spectra at two temperatures: 77 and 3007^ are demonstrated 
in the Figure 3, and the PLE spectra at those temperatures are depicted in 
the Figure 4 of the paper [Hj. They both show the dramatic changes of the 
band blue shift when the temperature is increased. Similar studies have also 
been performed before by the other groups [3 |H] . They observed the similar 
behaviors. 

Another interesting investigation is related to the temperature effect on 
the spectral line narrowing (SLN) of the poly(2-methoxy-5-(2'-ethylhexyloxy)- 
1,4-phenylenevinylene) (Mi^i/ — PPV) spin-coated from either THE or CB 
solvents |11 . In the experiment done by Sheridan et al., the SLN is mea- 
sured together with the absorption and PL as shown in Fig. 1 of their paper. 
It is found that the same trend of the SLN blue shift is observed as that for 
the absorption and PL with an increasing temperature. They attribute this 
to the same reason of the electronic structure modification resulting from the 
variation of the temperature. 

2.2 Temperature effects on structure, dynamics, and 
folding/unfolding of biomolecules 

Biomolecules are complex systems, featuring a large molecular size, a hetero- 
geneity of atomic constitutes and a variety of conformations or configurations. 
Their energy landscape thereby exhibits multiple substates and multiple en- 
ergy barrier, and varies in size for the barrier heights [23 1221 12SI 123 12H1 • The 
temperature should have a strong effect on their structure and dynamics in- 
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eluding the folding or unfolding [17-57]. This effeet eould be either from the 
fluctuation of thermal motion of molecules or due to the redistribution of 
electronic charge as we will discuss in the next section. 

The experimental observation of the temperature impact on the micro- 
scopic structure of biological systems dates back to the very early days. One 
focus, for example, is on the measurement of the elastic properties of the hu- 
man red blood cell membrane as a function of temperature ^l^j. Another 
related study is about the influence over the thermal structural transition of 
the young or unfractionated red blood cells due to the involvement of protein 
spectrin which might modified the spectrin- membrane interaction fl^ IHO] . 
Most recently, the atomic force microscopy {AFM) has been used to detect 
the impact from the variation of the temperature on the spectrin protein 
unfolding force as well as on the bond rupture force for the biotin-avidin 
system HH 1201 IHl • 

The AFM is a surface imaging technique with an atomic-scale resolution 
capable of measuring any types of the forces as small as 10"^*^ A^. It combines 
the principle of scanning tunnelling microscopes {STM) and stylus profilome- 
ter, and therefore can probe the surfaces of both conducting and nonconduct- 
ing samples [2I1E21- The imaging on soft materials such as biomolecules with 
the AFM has been performed beginning in the 1980's jHSlEllES]- Recently, it 
has been applied to measure the adhesive forces and energies between the bi- 
otin and avidin pair also as we will show in the next subsection [HHlEZlEHlEni- 
Unlike other experimental techniques, the AFM features a high precision and 
sensitivity to probe the surface with a molecular resolution, and can be done 
in physiological environments. 
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In an AFM investigation of mechanical unfolding of titin protein, for ex- 
ample, the restoring forces all show a sawtooth like pattern with a definite 
periodicity. It reveals much information about the mechanism of unfold- 
ing processes |421 143j . The observed pattern, in addition to fit a worm-like 
chain model, has also been verified by the steered molecular dynamics or 
Monte carlo simulations 003]. Similar study has been extended to other 
systems UHl iH HH] . 

The same kind of experiments has also been performed by varying the 
temperature. In the experiment carried out by Spider and Discher et al. pi] , 
the spectrin protein is chosen for the AFM study at different temperatures. 
Thousands of tip-to-surface contacts are performed for a given temperature 
because of the statistical nature of the AFM measurement. The observed 
curve for the relation between the unfolding force and extension shows the 
similar sawtooth pattern for all temperatures. In addition, the tandem repeat 
unfolding events are more favored at lower temperature as demonstrated in 
the unfolding length histograms. Most striking is that the unfolding forces 
show a dramatically nonlinear decreasing relation as the temperature T ap- 
proaches the transition temperature T^. This is shown in the Figure 3B of 
paper pT| . 

Similar behaviors regarding the force-temperature dependence have also 
been observed via either AFM or optical tweezers for the forced overstretch- 
ing transition for the individual double-stranded DNA molecules [121 1201 1121 

Some other interesting experiments which illustrate the effect of temper- 
ature on the microscopic structure of biomolecules have also been performed 
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even though the detailed physical origins of the effect (from either the elec- 
trons or the molecules) have not been specified EHl EH ESI EHl EZl ■ In 
a circular dichroism {CD) spectra and high resolution NMR study, for in- 
stance, it shows that the secondary structure of the Alzheimer jS (12-28) 
peptide is temperature-dependent with an extended left-handed 3i helix in- 
terconverting with a flexible random coil conformation [52,. Another example 
is related to the analysis of the temperature-dependent interaction of the pro- 
tein SshlOb with DMA which influences the DNA topology [SHlEllEni- The 
study from the heteronuclear NMR and site-directed mutagenesis indicates 
that the SshlOb exists as a dimer: T form and C form. Their ratio is deter- 
mined by the Leu^^ — Pro^^ peptide bond of the SshlOb which is sensitive to 
the temperature. 

2.3 Temperature effects on inter molecular forces 

The study of the general issues related to the temperature effects on the 
microscopic structure has been most recently extended to the realm of in- 
termolecular forces. Since the usual intermolecular forces such as hydrogen 
bond, van der Waals force, ionic bond, and hydrophobic interaction are weak 
and typically of the order of 0.1 eV or 4-0 kT at the physiological temper- 
ature, the variation of temperature will thereby have a very strong influence 
over the strength of these forces. 

The first experimental investigation on the temperature-dependent inter- 
molecular forces is for the biotin-avidin system and by an AFM measure- 
ment jnSl- The biotin-avidin complex is a prototypical receptor and ligand 
system with the biotin binding strongly up to four avidin protein jSHl EHl EHl 
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I5H |n2] • They have an extremely high binding affinity, and therefore serves 
as a model system for various experimental investigations. In the experi- 
ment carried out in Beebe^s group 63|, the receptor avidin is attached to the 
AFM tip and linked to the agarose bead functionalized with the biotin. The 
temperature of the entire AFM apparatus is varied at a range from 286 to 
310K. In addition, the loading rate is kept very slow so that the thermal 
equilibrium for the biotin-avidin pairs is assumed. The forces expected to 
be determined is the rupture force Fi between the individual biotin-avidin 
pair which is defined as the maximum restoring force jHSl- In actual AFM 
experiment, however, the total adhesive force between the tip and substrate 
is measured. It is a sum of finite number n of the interactions between each 
biotin and avidin pair. To extract the individual and average bond rupture, 
a statistical method has been developed in Beebe's group EUl IMl IM] . 
They assume a Poisson distribution for the number n of the discrete rupture 
forces or linkages from multiple measurements, and have obtained the single 
force Fi at different temperatures. The result is shown in the Figure 3 of 
the paper jUS]- We see that the individual rupture force Fi for the biotin- 
avidin pair is decreased by about five-fold in strength when the temperature 
is increased from 286 to 310K. 

To interpret the observed temperature impact on the biotin-avidin forces, 
Peebe's group has performed a thermodynamic analysis 63j. Based on the 
simple models and arguments [HZIEH], they have come out an equation that 
connects the square of the single bond-rupture force Fi to the absolute tem- 
perature T as follows. 




(1) 
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where the khond is the force constant of the individual biotin-avidin pair, and 
the time r^j is the characteristic time needed to break n pairs of those forces. 
The E"^ is the energy required to remove the biotin from avidin's strongest 
binding site and the corresponding time is t^. More detailed on this analysis 
can be found in the paper pJj. The relation between the square of the force 
Fi and temperature T is also plotted as Figure 5 in that paper. Therefore, 
from the relation (1) and that Figure, the information about the stiffness of 
the ligand and receptor bond and the critical binding energy, etc. can be 
obtained. Obviously, what we need at present is a microscopic theory which 
can account for all these relations and properties. 

2.4 Temperature effects on absorption spectra of elec- 
trons in solvents 

The structure and dynamics of solute in solvent is one of the most impor- 
tant fields in chemistry since most of the chemical reactions occur in solution 
phases. In the meantime, it is also one of the most challenging fields in the- 
oretical chemistry with many unsettled issues. The variation of temperature 
in the measurement of absorption spectra of solvated electron in various sol- 
vents has proved to be a useful means for the understanding of the solvation 
processes [71-85]. 

There are several experimental techniques available for this type of stud- 
ies with the pulse radiolysis being the most commonly used one. There are 
also several research groups doing the similar experimental investigations 
and obtaining the consistent results relating to the temperature effects on 
the optical absorption spectra of solvated electron in solvents. In a recent 
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experiment carried out in Katsumura^s group, for example, the pulse radi- 
olysis technique is employed to study the optical absorption spectra of the 
solvated electron in the ethylene glycol at different temperatures from 290 
to 598K at a fixed pressure of 100 atm. In addition to the faster decay of 
absorptions, it is found that, their maximal positions shift to the red with 
the increasing temperature as shown in the Figures 1 to 3 of the paper |HT]. 
This is in contrast to the situation for the electronic spectra of the polymers. 
They also point out the need to quantify the change of the density in the 
experiment in order to really understand the observed results. 

The same type of experiment has been extended to the study of the optical 
absorption spectra for Ag^ and Ag2 in water by varying the temperature, 
and similar results have been obtained [SUj . 

3 Theoretical Development 

Having presented four different types of experiments above, we have seen 
that the temperature effect on the microscopic structure of molecules is a 
very interesting and sophisticated field. Many more needs to be understood 
and probed. Even though the experimental study has been for a long time, 
very limited number of the related theoretical work is available, especially 
at the first-principle level. In other words, the quantum chemistry at finite 
temperature is not a well-defined or well-established field [U E] . 

It is true that the influence of temperature on the microscopic structure is 
a complicated phenomenon. There exists different functioning mechanisms. 
One consideration is that the variation of the temperature, according to 
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the Fermi-Dirac statistics, will change the thermal probability distribution 
of single-particle states for a free electron gas. It is expected that similar 
situation should occur for an interacting electron system, and therefore its 
microscopic structure will be correspondingly altered. Another consideration 
is that, for molecules or solids, the thermal excitation will cause the change of 
the time scales for the molecular motions. This will most likely bring about 
the transitions of the electronic states, and therefore lead to the breakdown 
of the Born-Oppenheimer approximation. Electron-phonon interaction is a 
fundamental topic in solid state physics and its temperature dependence is 
well-known. As a result, the temperature could change the strength of the 
coupling between the electronic and molecular motions. Nevertheless, we 
tackle the issues in a simpler way. We treat only an interacting identical 
fermion system, or neglect the coupling of the electronic motion with those 
of the nucleus in molecules or solids. We expect that some sort of the general 
conclusions will come out from this study. As a matter of fact, this is also the 
approach usually adopted in a non-adiabatic molecular dynamics, in which 
purely solving the eigenequation for the electrons will provide the reference 
states for the investigation of the coupling motions between the electrons and 
the nucleus of the molecules. 

In the following, we will present a self-consistent equation within the 
framework of equilibrium statistical mechanics which decides the molecular 
orbitals at a given temperature. 
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3.1 Hierarchy Bloch equations for reduced density op- 
erators in canonical ensemble 

We consider an identical and interacting A^-particle system. In a canonical 
ensemble, its A^th-order density operator takes the form 

= exp{-pH^), (2) 

and satisfies the Bloch equation [HUl E] 

- ^^"^ = H^D"", (3) 

where 

N N 

HN = J2hit) + Y.g{t,j), (4) 

i=l i<j 

is the Hamiltonian for the particle system composed of the one-particle 
operator h and two-body operator g. The (3 is the inverse of the product of 
Boltzmann constant kB and absolute temperature T. 

Since the Hamiltonian (4) can be written as a reduced two-hodj operator 
form, the second-order reduced density operator suffices to describe its (> 
2) particle quantum states. A pth-order reduced density operator is generally 
defined by jHS ES] 

= (5) 

where is the contraction operator acting on an A^th-order tensor in the 
A^-particle Hilbert space V'^ . The trace of the gives the partition function, 

Tr{DP) = Z{P,V,N). (6) 

Rewrite the Hamiltonian in a form 

N p N N 

j=p+l i=l j=p+l «<j(«>P+l) 
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where 

m = tK^) + t9{^,J), (8) 

1=1 i<j 

and apply the contraction operator on both sides of the Eq. (3), we 
develop an equation that the pth-order density operator satisfies |1] 



.1=1 



-^D^ = HfD^+{N-p)Ll^,[h{p+l)U^^'] + {N-p)Ll^, 



It provides a law according to which the reduced density operators vary in 
terms of the change of temperature. 

3.2 Hierarchy Bloch equations for reduced density op- 
erators in grand canonical ensemble 

The above scheme for deducing the equations for the reduced operators can 
be readily extended to the case of a grand canonical ensemble, which is a 
more general one with a fluctuating particle number A^. In this ensemble, 
the density operator is deflned in the entire Fock space 

oo 
N=0 

and is written as the direct sum of the density operators Dc{N) associated 
with the iV-particle Hilhert space V^^, 



Y,g{t,p+l)D^^' 



N=Q 



(10) 



+ 



(9) 



where 



Dg{N) = exp[-/3(if-^iV)], 
= exp(-/5^). 



(11) 



15 



and 



H = H-fiN, (12) 

is called the grand Hamiltonian on . The form of the Hamiltonian H has 
been given by Eq. (4) and the /i is the chemical potential. The corresponding 
pth-order reduced density operator is therefore defined as 



N=p 



P 



(13) 



with the trace given by 



Tr{D: 



N 
P 



and 



(14) 



(15) 



The E{P,fi, V) is the grand partition function. 

In a similar manner, we can also derive the hierarchy equations that the 
reduced density operators in the grand canonical ensemble obey j3] 



= H!D^+{p+l)L;^,[hip+l)D^+']+{p + l)Ll^, 



li=l 



J29{hP + 



+ 



p + 2 



TP 



gip+l,p + 2)DP+^ 



where 



and 



p 

i=l 



i<j 



(17) 



(18) 



h{i) = h{i) — /i. 

It gives us a law with which the reduced density operators in the grand 
canonical ensemble vary in temperature. 
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3.3 Orbital approximation and Hartree-Fock equation 
at finite temperature 

The Eqs. (9) and (16) define a set of hierarchy equations tliat establisli tlie 
relation among tfie reduced density operators D^, D^'^^, and Tfiey 
can be solved either in an exact scheme or by an approximate method. The 
previous study of N electrons with an independent particle approximation to 
the Schrddinger equation for their pure states has lead to the Hartree — Fock 
equation for the molecular orbitals [94-100]. We thereby expect that the same 
approximate scheme to the reduced Block equations (9) or (16), which hold 
for more general mixed states, will yield more generic eigenequations than 
the usual Hartree — Fock equation for the molecular orbitals. 

We consider the case of a grand canonical ensemble. When p = 1, Eq. 
(16) reads 



-^D' = H,D'+^^^^D'-^^D^hD'+2Ll [g{l,2)D^]+3Ll [g{2,3)D 

(19) 

Under the orbital approximation, the above second-order and third-order 
reduced density operators for the electrons can be written as 

D^^D^ AD^ A D^/{D^f (20) 

and 

L>2 = L»i A L>VL>°. (21) 

These are the special situations for the statement that a pth-order reduced 
density matrix can be expressed as a p-fold Grassmann product of its first- 
order reduced density matrices. With these approximations, the last two 
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terms of Eq. (19) can be evaluated in a straightforward way as follows 

2L\\g{l,2)D^\ = {J-K)D\ (22) 



and 



where 



and 



3L^ 



5(2,3)D= 



Tr{gD 



1 



J^Tr, 



'g-D\2-2)\ /D\ 



(23) 
(24) 



g-{2,^)-D\2-2)\/D\ (25) 

are called the Coulomb and exchange operators, respectively. With the (2, 3) 
being the exchange between the particle 2 and 3, the action of the K on the 
reduced density operator is 



/s:-Di(3;3) = Tr2 (2,3) ■Di(2;2)l /D°-Di(3;3) 



g ■ D\'i-2) ■ D\2-M /D 



(26) 



Substitution of Eqs. (22) and (23) into Eq. (19) yields the Block equation 
for the first-order reduced density matrix of interacting electrons under 
orbital approximation, 



•9 ^1 (TrhD^ TrD^ 



\ - ^D\F - ^)D\ (27) 



where 

F^h + J-K, (28) 

is called the Fock operator at finite temperature. Redefine the normalized 
first-order reduced density operator 

pi = DV^o, (29) 



we can simplify above equation into 

--^p'^{F-n)p'-p\F-n)p\ (30) 

Furthermore, from Eq. (30) and its conjugate, we get 

Fp'-p'F^O, (31) 

which means that the Fock operator F and the first-order reduced density 
matrix commute. They are also Hermitian, and therefore have common 
eigenvectors These vectors are determined by the following eigen 

equation for the Fock operator, 

=e#,). (32) 
The first-order reduced density operator is correspondingly expressed as 

= X]^(/?,)",ej)|0i)(0i|, (33) 

i 

where a;(/9,/i, e^) is the thermal probability that the orbital is found to be 
in the state at finite temperature T. Substituting Eq. (33) into Eq. 

(30), we can obtain the equation this thermal probability a;(/3, p, ei) satisfies, 
d 

Its solution has the same usual form of the Fermi — Dirac statistics for the 
free electron gas as follows, 

uj{(3,p,ei) = -^^ (35) 

with the energy levels {ej} determined by Eq. (32). 
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4 Discussion, Summary and Outlook 



In this paper, we have presented a description of both experimental and the- 
oretical developments related to the temperature impact on the microscopic 
structure and processes for the molecules. 

In the theoretical part of this paper, wc have depicted the sets of hier- 
archy Block equations for the reduced statistical density operators in both 
a canonical and a grand canonical ensembles for the identical fermion sys- 
tem with two-body interaction. We have solved the equations in the latter 
case under a single-orbital approximation and obtained an eigen-equation 
for the single-particle states. It is the extension of usual commonly used 
Hartree — Fock equation at the absolute zero temperature to the situation 
at any finite temperature. The average occupation number formula for each 
single-particle state is also obtained, which has the same analytical form as 
that for the free electron gas with the single-particle state energy determined 
by the Hartree — Fock equation at finite temperature (32). 

From Eqs. (24), (25) and (28), we see that the Coulomb operator J, 
the exchange operator K, and therefore the Fock operator F are both the 
coherent and the incoherent superpositions of single-particle states. They 
are all temperature-dependent through an incoherent superposition factor, 
the Fermi — Dirac distribution, uj{/3, iJ,,ei). Therefore, the mean force or the 
force field, and the corresponding microscopic structure are temperature- 
dependent. 

We have expounded the physics relating to the temperature effects on 
the electronic structure or processes of the molecules. This is very critical for 
our understanding and studying the temperature infiuence over the molecular 
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structure. From this analysis, for example, we could see that the temperature 
should have a stronger effect on the molecular transition states and therefore 
their chemical reactivity. Starting from Eqs. (9) or (16), it will be a very 
significant work to establish a corresponding multireference theory for the 
molecular orbitals at finite temperature jHEl HHH 11021 • 

On the experimental sides, we have exposed four major fields of investi- 
gations of chemical or biomolecular interests, which show the temperature 
impact on their structures, spectra, or bond rupture forces. 

The complete determination of the geometric structure and electronic 
spectra of the polymeric molecules is a very difficult task. As has been 
stated in papers fT^ I16j . there are many different components contributing 
to the change of the spectra. At present, we focus on the study of the effect 
from the temperature. We have demonstrated that it can alter both the 
shapes and positions of the absorption and other spectra for the PPV and 
its derivatives. As has been analyzed, the increase of the temperature will 
cause the excitation of the vibrational, rotational and liberal motions which 
might also lead to the electronic transition. The Huang — Rhys parameter 
has been introduced to describe the strength of the coupling between the 
electronic ground- and excited-state geometries. Furthermore, it has been 
found that this factor is an increasing function of the temperature [3 |H1 
p. Obviously, a more detailed study of the electronic structure, excitation 
and spectroscopic signature at the first-principle level which includes the 
temperature-dependent force field is expected. 

The temperature has proved to be a big player in both experimental and 
theoretical study of the structure and dynamics of biomolecules including 
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their folding or unfolding. At a first glance, the energy gap between the 
HOMO and LUMO for the biomolecules should be small or comparable to 
the Boltzmann thermal energy fc^T because of their very large molecular 
size. Therefore a change of the temperature should have a strong influence 
over their electronic states, and consequently, the energy landscape and the 
related dynamics including the folding or unfolding, etc. The experimental 
investigation with the AFM and other techniques of the temperature effect on 
the shift of their unfolding pathways might have verified this sort of thermal 
deformation of potential energy landscape 1201 1^ • This is in contrast to 
the tilt and deformation of energy landscape of biomolecules including their 
transition states resulting from the applied mechanical forces |69j . 

Unfolding proteins by temperature is not just one of the classical exper- 
imental techniques for the study of the structure, dynamics and energetics 
of the biomolecules. It has also been employed, for example, in the molecu- 
lar dynamic simulation to study the structure of transition states of CI2 in 
water at two different temperatures: 298 K and 498 K [22] • The later high 
temperature is required in order to destabilize the native state for monitoring 
the unfolding as done in the real experiments. In another recent molecular 
dynamics simulation Karplus's group has compared the temperature- 
induced unfolding with the force-stretching unfolding for two /^-sandwich 
proteins and two a-helical proteins. They have found that there are the sig- 
nificant differences in the unfolding pathways for two approaches. Neverthe- 
less, in order to get more reliable results, temperature-dependent force fields 
need to be developed and included in the molecular dynamics simulations. 
This is also the case in the theoretical investigation of protein folding since 
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an accurate simulation of protein folding pathways requires better stochastic 
or temperature-dependent potentials which have become the bottleneck in 
structure prediction [211 123 123 123 12H1 ■ From structural points of view, the 
variation of temperature leads to the change of the mean force or the energy 
landscape, and therefore provides a vast variety of possibilities, for instance, 
in the protein design and engineering. 

The intermolecular forces are ubiquitous in nature. They are also ex- 
tremely important for the biological systems and for the existence of life. 
The intermolecular forces have the specificity which are responsible for the 
molecular recognition between receptor and ligand, antibody and antigen, 
and complementary strands of DNA, and therefore for the regulation of com- 
plex organization of life |lU3j . For these reasons, the experiment carried out 
in Beebe's group has an immediate significance. It has demonstrated that 
the temperature can be an important factor for changing the specificity of 
the intermolecular forces and therefore the function of life jU^ . Nevertheless, 
how the charge redistribution occurs due to the variation of the tempera- 
ture has not been interpreted, and a microscopic theory for quantifying the 
temperature influence on the intermolecular forces is still lacking. Since the 
delicate study of the intermolecular forces provides the insight into complex 
mechanisms of ligand-receptor binding and unbinding processes or pathways, 
a paramount future research is to establish the links between the intermolec- 
ular forces and the temperature within the quantum many-body theory. 

The theoretical study of the temperature effects on the optical absorp- 
tion spectra of solvated electron is still in very early stage and few published 
works are available jH2l EHl Ell ESI- One of the earliest studies by Jortner 
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used a cavity model to simulate the solvated electron where the electron is 
confined to the cavity surrounded by the dielectric continuum solvent jH2]- 
However, his study is not of fully microscopic in nature since he assumed 
a temperature dependence of phenomenological dielectric constants which 
were also obtained from the available experimental data. In addition, the 
model used is too simplified and, for instance, it neglects the intrinsic struc- 
ture of the solvent molecules. There are a few recent investigations on the 
temperature effects on the absorption spectra of solvated electrons. They all 
cannot catch the full features of the experimental observations. One reason 
is that the physical nature for the process is not totally understood which 
might leads to incorrect models used for the simulation. The other is to 
utilize the crude models which might have omitted some important physical 
effects. For example, in an analysis by Brodsky and Tsarevsky [23] , they have 
concluded a temperature-dependence relation for the spectra which is, how- 
ever, in contradiction with the experimental findings at high temperature. 
The quantum path-integral molecular dynamics simulation cannot produce 
those temperature-dependence relations observed in the experiments [H5- In 
a recent quantum-classical molecular-dynamics study by Nicolas et al, even 
though the temperature-dependent features of optical absorption spectra for 
the solvated electron in water have been recovered |H3], however, they claim 
that the red shifts of absorption spectra with increasing temperature observed 
in both experiments and calculation are due to the density effect instead of 
temperature. This might cast the doubt of usefulness of our current theoret- 
ical work in this area. However, after examining their work, we observe that 
they actually have not included any temperature effect on the electron in 
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their theoretical model. This effect might be either from the Fermi — Dirac 
distribution for individual electrons or due to the electronic excitation caused 
by the thermal excitation of the solvent, as we have discussed in paper 
Obviously, much finer theoretical work or more experimental investigations 
in this area are expected to resolve this dispute. 

In addition to the systems discussed above, there are other types which 
also show the temperature impact on their microscopic structures. Either 
theoretical or experimental work have been done or are in progress. Examples 
include the study of the temperature dependence of the Coulumb gap and 
the density of states for the Coulomb glass, the experimental investigation of 
the temperature effects on the band-edge transition of ZnCdBeSe, and the 
theoretical description of the influence from the temperature on the polaron 
band narrowing in the oligo-acene crystals |lU4t I1U51 llU6j . 

To sum up, the quantum chemistry at finite temperature is a new and 
exciting field. With the combination of the techniques from current quantum 
chemistry with those developed in statistical or solid state physics, it will 
provide us with a myriad number of opportunities for the exploration. 
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